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I visually explore the features of geodesic orbits in arbitrary stationary axisymmetric vacuum 
(SAV) spacetimes that are constructed from a complex Ernst potential. Some of the geometric 
features of integrable and chaotic orbits are highlighted. The geodesic problem for these SAV 
spacetimes is rewritten as a two degree of freedom problem and the connection between current 
ideas in dynamical systems and the study of two manifolds sought. The relationship between the 
Hamilton- Jacobi equations, canonical transformations, constants of motion and Killing tensors are 
commented on. Wherever possible I illustrate the concepts by means of examples from general 
relativity. This investigation is designed to build the readers' intuition about how integrability 
arises, and to summarize some of the known facts about two degree of freedom systems. Evidence 
is given, in the form of orbit-crossing structure, that geodesies in SAV spacetimes might admit, a 
fourth constant of motion that is quartic in momentum (by contrast with Kerr spacetime, where 
Carter's fourth constant is quadratic). 
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I. INTRODUCTION 

The Carter constant associated with a Kerr spacetime 
plays a crucial role in current LIGO (Laser Interferom- 
eter Gravitational Wave Observatory) and LISA (Laser 
Interferometer Space Antenna) extreme and intermedi- 
ate mass ratio inspiral (EMRI/ IMRI) waveform calcula- 
tions. Its generalization to arbitrary stationary axisym- 
metric vacuum spacetimes (SAV), could lead to an al- 
gorithm to map spacetimes around a compact object de- 
scribed by a general set of multipole moments; see [l[ and 
Paper I of this series . Such an algorithm may pro- 
vide a method of determining the nature of compact ob- 
jects by asymptotically observing gravitational and elec- 
tromagnetic radiation from an EMRI or IMRI inspiral. 

The geodesic problem in Kerr spacetime is completely 
solved by the specification of four isolating integrals or 
constants of motion, namely rest mass, energy, axial an- 
gular momentum and the Carter constant (//, E, Lz,Q). 
The generalization of the first three constants of motion, 
namely (/i, E, L^) to SAV spacetimes is trivial. These 
constants result from the absence of an explicit depen- 
dence of the Lagrangian on proper time, coordinate time 
and the axial angular coordinate respectively (r, t, (/>). 
The meaning of the fourth constant Q, first discovered 
by Brandon Carter by separation of the Hamilton- Jacobi 
equations (HJE) 0, Q is a little more obscure. It is this 
fourth constant that allows the reduction of the geodesic 
equations to first-order quadratures, and the complete 
solution of the geodesic problem. The geometric inter- 
pretation of this fourth constant and the conditions for 
its existence are explored in this and subsequent papers 
in this series [si, . 

The present paper visually characterizes the geodesies 
in some of the axisymmetric spacetimes that are gener- 
ated from an Ernst potential. In particular, the Manko- 
Novikov 0, H, [i] and Zipoy-Voorhees metrics [13, [H 
are considered. The geodesic problem is formulated as 
a two degree of freedom problem (2-DOF) in dynami- 



cal systems. Ideas from the field of integrable systems 
are collated and introduced by means of a series of vi- 
sual examples. For historic purposes, the role of the 
Hamilton- Jacobi equation (HJE) is put in context. Pos- 
sible tests for integrability are addressed. The concepts 
of phase and energy space are introduced and illustrated 
by means of an example. The role and possible forms 
of the additional invariant are explored and a geometric 
interpretation of Killing tensors given. 

Finally, some of the frustrations and computational 
difficulties when dealing with 2-DOF Hamiltonians are 
mentioned, and the implications of the numerical exper- 
iments in SAV spacetimes for the existence of a general- 
ized Carter constant are described. 



II. EQUATIONS OF MOTION 

I begin with the general SAV spacetime line element 
of the form 



o2V 



{dt 
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where i/', 7, ^ Eind R are functions of p and z and k 
is a real constant. The vacuum field equations relate 
these functions to solutions of the Ernst equation for the 
complex potential £, 
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where V 
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^pp ~^ p^p 
dot is the flat-space inner product 



V — {dp,dz), and the 
In particular, the 



function e^^ = 5R(£) denotes the real part of the poten- 
tial. The functions 7 and uj can be obtained by means of 
line integrals of the potential once it is known. A gauge 
freedom in the form of the harmonic function R obeying 
Rzz + Rpp = exists in this metric. Often this freedom is 
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used to set R= p: however, for the sake of later compar- 
ison with solution generation techniques, I shall retain 
this generality. 

The Hamiltonian associated with geodesies of this met- 
ric is 



1 



(3) 



where following the notation of Goldstein 

q — {p, z, (f), t) are the generalized coordinates and 

p = (jjp, Pz, Pci>, Pt) are the conjugate momenta. 

In order to write the equations of motion in compact 
form I make use of the Poisson brackets. The Poisson 
bracket of two functions g and h with respect to the 
canonical variables is defined as 



dg dh 
dqk dpk 



dg dh 

dpk dqk 



(4) 




The geodesic equations can now be expressed in first or- 
der form using Hamilton's equations, namely 



Pti = [Pf^,'H], 



(5) 



where the dot • indicates the total derivative with respect 
to proper time r. 

Using this notation it is immediately obvious that the 
absence of any explicit metric dependence on t and 
results in Pt = p^ = 0. By setting these quantities equal 
to the standard constants pt — —E and p^p = Lz, the 
study of geodesic motion in four-dimensional spacetime 
is reduced to the study of a 2-DOF dynamical system 
with an effective potential. The reduced Hamiltonian 
can be expressed as: 



H{p,z,pp,pz) 



Pz 



G 



(6) 



where the two potentials V and G have been introduced 
to simplify notation, and are defined as 



V{p,z) = ee''^~'^, 
G{E,Lz,p,z) 
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with A, B indicating the components i, 0, and let i,j 
range over p, z. The Hamiltonian constant H = —1/2 p? 
fixes the sum of the squares of the conjugate momenta to 



pI +pl^iG~ p^)V = J(p, z, E, L, m"). 



and the equations of motion become 



Pi 

q. = 



Pi 



2V ■ 



(9) 
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Upon introducing the non afhne parameter A such that 
VdX — dr, and letting ' indicate differentiation with re- 
spect to A, the equations further simplify to 



Pi 



Pi 



(11) 



FIG. 1: Constant J potential surfaces for a geodesic in the 
Schwarzschild metric with E = 0.95, Lz = 3, /i = 1. Contour 
Spacing 0.5. Explicit form of the metric functions can be 
found in Appendix [Cl 



The problem of finding a generalized Carter constant 
in SAV spacetimes can be expressed most generally as 
the hunt for a function Q{p, z,pp,pz) distinct from the 
Hamiltonian H{p, z,pp,pz), that remains constant along 
an orbit of the 2-dimensional Hamiltonian H , i.e. a func- 
tion such that [Q,H] — 0. Alternatively this can be 
stated as, the study of the geodesies of the two manifold 
with metric gj and associated "Jacobi" Hamiltonian Hj, 



dJij — J5ij, 



Hj = 



(pI+pD 

2J 



(12) 



A more general and rigorous treatment of these ideas is 
given in |13l |. 

In general, for a generic 2-dimensional Hamiltonian 
Hg, no such integral of motion Qg exists and the Hamilto- 
nian is chaotic. In most textbooks on dynamical systems 
completely integrable systems are given but brief men- 
tion [13, liljIlBl- Explicit examples are rare. A thorough 
review summarizing most of the known examples can be 
found at p^ . 



III. A FEW WORDS ABOUT THE HAMILTON 
JACOBI EQUATION 

Carter's original derivation [3)] of the fourth invariant 
for the Kerr metric was performed by means of separation 
of the Hamilton Jacobi equations (HJE) for the Jacobi 
function S, 



1 p,d^d^ 
2^ ^q^^^q•'' 



(13) 
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The Jacobi function generates the canonical transforma- 
tion to action angle variables [T^ fPage 449). Once S is 
known, the problem of finding the full set of constants of 
motion is solved. One method of solution of and the 
only one so far used in practice, is by means of separation 
of variables. 

Lack of separation of variables, however does, not nec- 
essarily imply that the system is not integrable or the 
absence of Q. An example that does not have its origin 
in the vacuum field equations is the Fokas-Lagerstrom 
Hamiltonian 



which admits an orbital invariant 



-2/3 



Q = ipl 



(14) 



-2/3 

(15) 



No separation of variables for the HJE associated with 
has ever been found. 
In the case of the SAV spacetimes, all metrics admit- 
tinga second-rank Killing tensor have separable HJE's 
[1, 13| in some coordinate system. This feature will be 
considered in greater detail in Paper III of this series Q 
where I will catalog the coordinate systems where this 
occurs. 
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Many of the difRculties in carrying out a test are sum- 
marized in [IQ] and some will be demonstrated later in 
this paper. Partial tests, such as the Painleve test, can 
be carried out and all dynamical systems that pass this 
test have been found to be integrable. Failure to pass 
the Painleve test however does not imply that a partic- 
ular Hamiltonian will fail the test in another coordinate 
system. 

The Zipoy-Voorhees Metric |l3, Hii fails the Painleve 
test in the same manner that the Fokas-Lagerstrom 
Hamiltonian expressed in the form (|14p does. This fact is 
inconclusive since after a number of difficult transforma- 
tions a formulation of the Fokas-Lagerstrom Hamiltonian 
was found that passes the Painleve test [l^ . 

Possibly the strongest indication that integrability fails 
is the finding that numerical integration yields a Poincare 
map without closed curves. A Poincare map that displays 
closed curves is indicative that the Hamiltonian may be 
integrable, but it does not provide proof of the existence 
of an additional constant of motion. 

If a Hamiltonian H is close to an integrable Hamil- 
tonian Hq with invariant Qpi it is always possible, fol- 
lowing a perturbative scheme developed by Deprit [l9|, 
to compute an approximate invariant Q associated with 
H, and thus to produce an approximate Poincare map 
with closed curves. A perturbative invariant so con- 
structed, however, may not give an accurate rendition 
of the phase space of the perturbed Hamiltonian iJ in a 
strongly chaotic regime. The classical example where this 
is clearly illustrated is the Henon-Heiles problem [13] . A 
perturbative analysis is not sufficient to prove or disprove 
integrability. 

Attempting to use the Deprit scheme of canonical per- 
turbation theory to construct invariants for SAV space- 
times is prohibitively expensive computationally, and is 
not feasible if a solution for all SAV spacetimes is sought. 

Indications of integrability can also be gleaned by ob- 
serving the structure of the orbits in configuration space. 
This will be discussed in greater detail in subsequent 
paragraphs. 

In Paper IV of this series @ I will propose a test to see 
if SAV spacetimes admit invariants that are polynomial 
in momenta. 
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ORBITS, PHASE SPACE, AND ENERGY 
SPACE 



FIG. 2; Configuration (p, z) space depiction of a geodesic 
orbit for the potential J shown in Figure [T] (Schwarzschild 
metric with E = 0.95, = 3, /i = 1) 



IV. TESTS FOR INTEGRABILITY 

To my knowledge, there exist no conclusive algebraic 
tests for integrability for a given 2-DOF Hamiltonian H. 



Integrable systems have a surprisingly simple struc- 
ture [14| . If expressed in terms of action- angle variables, 
the orbits are found to trace out tori in the four dimen- 
sional phase space {p, z,pp,pz). Although four dimen- 
sions are difficult to visualize, it is possible to see what 
these orbits look like in the three dimensional energy 
space by introducing a momentum phase angle 6 such 
that 



Pp 



— Vj cos 9, 



P^ = VJ sine. (16) 
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Doing so explicitly imposes, the Hamiltonian con- 
straint ([5]), so the orbit can be visualized in {9, p, z) 
energy space, as depicted in Figure [3l Note that if the 
light blue lines in this figure are "squashed" , i.e. pro- 
jected down onto the (p, z) plane, a rotated version of 
Figure [2] is obtained. 




FIG. 3: Three dimensional energy space depiction of the torus 
of the orbit displayed in Figure (2] The orbit is depicted in 
cyan lines. The surfaces depicted correspond to the surfaces in 
phase space at which the orbit in configuration space reaches 
an extremum in curvature. 

One method of characterizing a curve in the (p, z) 
plane (Figure H]) that is independent of the parameter- 
ization of the curve is to compute its curvature k. (The 
curvature of a curve is a measure of how rapidly the curve 
is moving away from its tangent line.) 

The curvature of a curve parameterized by (p(t), z{t)) 
can be expressed as 



zp — pz 



(p2+i2) 



3/2 ■ 



(17) 



Using the momentum phase angle 0, this simplifies to 

K = i (sing dp[\nJ)-coBe d^{hiJ)). (18) 

There are special points along a curve at which the cur- 
vature stops changing, or reaches an extremum, namely 
when k — Q. These points are indicated by means of a 
dark blue line in Figure [2] The surfaces in energy space 
on which these extrema occur can be computed for any J 
(See appendix E]) . The brown surfaces shown in Figure 
13] are the extreme curvature surfaces. 

The rate of change of the phase angle along a particular 
geodesic can be calculated and expressed compactly in 
terms of the curvature as 



or 



(19) 



If the geodesies are integrable, as is the case in Fig- 
ures [5] and [21 the orbits sweeps out a surface in the en- 
ergy space (Figure [3]) . The locus of points at which the 
orbit reaches a point of extreme curvature forms a curve 
in configuration space (Figure [2]) . The points of contact 
with the J — Q contour are unique, and are determined 
by the constant Q. Furthermore, if a Poincare map is 
drawn it is constituted out of closed curves. 

If the geodesic problem is not integrable (no Q exists) 
and the orbit is strongly chaotic, it will wander all over 
energy space. If the geodesic problem is not integrable 
and the Hamiltonian H is close to an integrable Hamil- 
tonian, Hq, the orbit can be confined to a small volume 
in phase space, and appear integrable. If both H and i/o 
meet certain criteria |14| . the manner in which the sur- 
faces in energy space, or tori in configuration space, that 
exist for are destroyed are quantified by the KAM 
theorem 



VI. METHODS FOR CONSTRUCTING AN 
ADDITIONAL INVARIANT Q 

Suppose now that you suspect that an additional in- 
variant Q exists, due to a numerical exploration that 
yielded a Poincare map with closed curves or due to the 
fact that the points of extreme curvature all lie on a curve 
or due to your Hamiltonian passing the Painleve test; and 
suppose you would like to construct an explicit expres- 
sion for Q. The rules of engagement to date appear to be 
simple: you guess its form and hope that you are right. 
One method of guessing is to postulate that the invari- 
ant Q is polynomial in momenta p. This is equivalent to 
guessing that you have a Killing tensor on your two man- 
ifold. This is by no means the only form an invariant can 
take; however most of the known examples of integrable 
two dimensional Hamiltonians have polynomial Q's jT6| . 
(For a further discussion on the generality of this form of 
guessing see [2l|.) 

The task of finding an additional invariant Q poly- 
nomial in momenta for a 2-DOF problem in dynamical 
systems has a long history. In some avenues of litera- 
ture it is know as Whittaker's problem [li,|2l|. Hall [2l| 
provides a very complete and readable reference and the 
analysis adopted here is guided largely by his treatment 
of the problem. This analysis will become particularly 
useful in a Paper IV Q in this series when I return to 
a four dimensional representation of the geodesic prob- 
lem, and attempt to understand the coupling between 
the Weyl tensor of SAV spacetimes and the possible ex- 
istence of Killing tensors. This method of analyzing the 
additional invariant appears to identify the most impor- 
tant quantities that should be considered, and provides a 
geometric picture of what they are. Some of the difficul- 
ties in checking why a given Hamiltonian is integrable are 
also illustrated. Furthermore, this approach highlights 
other properties an integrable orbit has. The problem 
was also considered by [23| and in different notation 
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by I24|,l25|,l26|,l27|. 

Before I begin the analysis it is useful to introduce the 
complex variable C — l/2(/3 + iz). Let ( = l/2{p — iz) 
denote it's complex conjugate. In terms of this complex 
variable the orbital curvature in the (p, z) plane can be 
expressed as 



(20) 



and derivatives along the geodesic parameterized by A 
become 



d\ = Ppdp+Pzdz 



(21) 



Let us exploit the phase angle introduced in Eq. (|16p to 
express our additional invariant Q, which is assumed to 
be a general A'^th order polynomial in the momenta, pp 
and pz , as [cf. Eq. p6)) ] 



1 ^ 



(22) 



where the Q„ are complex valued functions of the con- 
figuration space variables, Qn = ^^(CtOi is a positive 
integer and Q-n — Qn- In effect, we are building up 
a Fourier series representation of the surface the orbit 
sweeps out in energy space. Figure |31 

The condition that Q is invariant along the orbit, in 
other words that Q' = 0, results in differential equations 
for the functions Qn- Explicitly, computing 9a (1221) and 
making use of 6' = — ka/J yields, 

n ^ 
n ^ 

+ ^ ^""/'5^(QnJ"/')e'("-^)''. (23) 

If this expression is to hold for all 9, the coefficients of 
gifce Yonst vanish for all -(iV + 1) < A: < (A^ + 1), which 
translates into the conditions. 



Qr. 



= 0, 



\ jn/2 



Qn-1 

'c I 

J 2 



for n = N,N-1, 

ioT 0<n< N ~l, 
(24) 



where Qo is real and ^{d^{Qi^/j)) = 0. For negative 
n values we get the complex conjugates of the above ex- 
pressions. These equations can be identified directly with 
the Killing equations for a two-manifold. The correspon- 
dence is shown in detail in Appendix [B) The lessons 
learned here will be exploited in that setting in my Pa- 
per IV 0. 



Just as in a Fourier decomposition, the equations for 
odd and even n decouple. Furthermore, the first condi- 
tion of (|24p implies that 



■N 



(25) 



is an analytic function of ^ and indicates an inherent 
"gauge" freedom in the Killing equations. It is this free- 
dom that makes the identification of integrable 2-DOF 
Hamiltonians so difficult. The integrability conditions 
that the functions Qn exist result in conditions on the 
conformal factor J. As a result, you can write down the 
differential equations the conformal factor J must obey if 
it is to admit a Killing tensor in some coordinate system. 
However if you are given a sample Hamiltonian to check 
for integrability, you have no idea what transformation 
leads to the coordinate system where we can conduct the 
check. An additional difficulty is that the conditions on 
the conformal factor for N > 2 are highly nonlinear. 

In the case of SAV spacetimes that admits a second- 
rank Killing tensor (Carter spacetimes) it is possible to 
exploit the coordinate freedom to our advantage by cou- 
pling it to the gauge freedom in the metric (The R func- 
tion in equation ([T])). A derivation of Carter spacetimes 
using this method is given in Paper III of this series . 

The case where N = 1, i.e. where the invariant Q is 
linear in the momenta, corresponds to a two-metric with 
conformal factor J that admits a Killing vector. This 
implies that it is a manifold of constant curvature (of the 
two manifold, not the orbit). The curvature of the two 
manifold is given by 



K 



1 



%(lnJ) 



(26) 



The three possibilities include flat space {K = 0), the 
two sphere {K > 0) and the Lobachevskii plane {K < 0) 
which can be visualized as the surface of a bugle. 

The problem of an invariant quadratic in the momenta 
{N — 2), on a two manifold was solved by Koenigs [1^ 
in 1889, who distinguished four types that are closely 
related to the four separable coordinate systems found 
by Carter (a derivation is given in Paper III [|) and 
to the super-integrable systems studied by Kalnins et al 
[H, [2^. The algebraic properties of two manifolds of 
this type are classified in [2J, [2^ . Koenigs provided a 
very accurate geometric description of what the second 
quadratic invariant actually represents. This geometric 
picture was revisited and generalized by Moser [28| , and 
clearly illustrated by Knorrer [2^ in his study of geodesic 
fiow on an ellipsoid. The second invariant corresponds 
to the Hamiltonian constant on a two manifold distinct 
from the first and there exists a very simple geometric 
construction mapping the geodesies on the one manifold 
to the next. 

For the = 4 case very few examples are known 
[1^ . [30| . It is my thesis that a large class, of these two 
manifolds are generated by SAV spacetimes. I further 
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suggest that the quartic structure is very closely related 
to the algebraic structure of the Weyl tensor, and that 
solution generation techniques for two manifolds already 
exist in the form of the solution generation techniques for 
SAV spacetimes. A test to see whether the SAV space- 
times admit a fourth order Killing tensor is proposed in 
Paper IV of this series [6| . Numerical evidence that indi- 
cates that SAV spacetimes might generate two manifolds 
with fourth order Killing tensor is given in Section IVIIl 



VII. NUMERICAL EXPERIMENTS AND 
ORBITAL STRUCTURE 

The existence of an invariant Q of the form of Eq. ([22)) . 
equivalently a Killing tensor, has direct implications for 
the orbital appearance of a geodesic. Consider a specific 
point in configuration space, for example a point in Fig- 
ure [H and consider the possible tangent directions that 
the geodesic could have leaving that point. The invari- 
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FIG. 4: An orbit in the Zipoy-Voorhees metric with 5 — 2. 
Orbital parameters are E — 0.95, L = 3, fi = 1. The black 
dots '•' and blue stars indicate extrema in orbital curvature. 

ant Q restricts the possible exiting tangent directions to 
the number of zeros of Eq. ([22|) . If = 4 and the orbit 
can be traversed both ways, as in the SAV case, the an- 
swer is that there are only 4 possible crossing directions 
(two crossing curves). This statement must hold for ev- 
ery point in configuration space. As a result, the orbits 
have a very ordered, cross-hatched appearance. Many 
numerically explored SAV spacetimes display this cross 
hatching pattern. 

One example for which this orbital structure is ob- 
served for all parameter values I explored, is the Zipoy- 
Voorhees Metric [13, EH- A special case of the Weyl class, 
this metric has the multipole structure of a finite rod. 
The metric functions are given explicitly in Appendix ID] 



It represents the one parameter (S) family of spacetimes 
that links flat space {5 — 0) to the Schwarzschild solution 
{S = 1). The orbital structure is displayed in Figure IH 
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FIG. 5: Metric components of the Manko-Novikov spacetime. 
Thick black contours indicate zero values, red contours or re- 
gions marked by > admit positive values and blue contours 
or < 0, negative. 

An example of possibly greater astrophysical applica- 
tion in the EMRI problem is the Manko-Novikov space- 
time ^, whose metric components are given in Ap- 
pendix [Cl The initial exploration into the orbits of this 
spacetime was performed by Gair et al. (sij and unusual 
orbital behavior was observed. 

Some of the properties of the Manko-Novikov space- 
time are sketched in Figure [10] and the functional form 
of the metric functions given in appendix [C] Figure \E\ 
characterizes the nature of the metric functions close to 
the horizon. 

In large regions of the parameter space of this space- 
time, for example region B of Figure [HI the geodesic or- 
bits display the characteristic fourth order crossing struc- 



1.5 2 2.5 3 3.5 4 

P 

FIG. 6: Two orbits in the Manko-Novikov spacetime. The metric functions for this spacetime are given at the end of AppendixICl 
and the following parameters were used a — 0.626789, 02 ~ 11.4708 and k — (1— a^)/(l+Q!^). The orbital parameters associated 
with the orbit in the figure are E = 9.5, L = —3 and jj. — 1. 



ture. There are however regions, discovered by Gair et 
al. [3l|, where integrability fails and the orbit is chaotic. 
One such region in which chaos occurs is displayed in Fig- 
ure [HI region A. The inset provides the contours of the 
potential function Ja for this region; it lies outside the 
ergoregions displayed in Figure [5] The Poincare maps 
for this orbit fail to display closed curves [3l| and the 
random orbital crossing structure implies that a Killing 
tensor on this manifold will never be found. I know of 
no similar example, in the literature of 2-DOF dynami- 
cal systems, where the orbits appear entirely integrable 
in one region and chaotic in another. 

The numerical experiments conducted by Gair et al. 
[3l| lead them to conclude that inspiralling orbits are un- 
likely to sample the "chaotic" region, so the possibility 
of observing such orbital behavior during a gravitational 
wave inspiral event is small, a conclusion with which 
I concur. However, conventional wisdom holds that if 
one observes the failure of integrability in some region 
of phase space it should preclude the construction of an 
invariant for the Hamiltonian in another. It would be an 
unfortunate and strange irony if "chaotic" behaviour in a 
region of phase space that is observationally inaccessible 
prevents us from obtaining an explicit expression for an 
invariant in the region of phase space from which observ- 
able gravitational radiation results. It is this quantity 
that will give us theoretical power in describing inspi- 
ralling orbits in an algorithm for mapping spacetime. 

Since the two regions A and B are disjoint they can 



be considered as two separate two manifolds J a and Jb , 
and following the analysis performed in this paper there 
is nothing that implies the chaos observed in region A 
precludes the existence of a Killing tensor on the two 
manifold Jb- To date the origin of the chaos in region 
A has not been carefully characterized. It is unclear 
whether the KAM theorem can be applied to this case, 
as the region A has no counterpart in the integrable Kerr 
spacetime, to which it reverts if the anomalous multipole 
moments are set to zero. In many ways the explanation 
of the orbital behavior in region A remains a very inter- 
esting puzzle. 



VIII. CONCLUSION 

This paper formulates the problem of finding the 
fourth invariant, more precisely, the generalization of 
Carter's constant to all SAV spacetimes, as a 2-DOF 
problem in dynamical systems. Equivalently stated, the 
problem can be formulated as the study of geodesic flow 
on a two manifold admitting a two metric with conformal 
factor J. A combination of the original metric functions 
in Eq. ([T|) and the constants that can be easily obtained 
from the metric symmetries determine J (equation ©). I 
summarize some related developments in dynamical sys- 
tems and in the study of two manifolds, which may help 
this problem and point out some of the difficulties faced. 
In particular, I emphasize the absence of a conclusive 
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algebraic check of whether a two-manifold is integrable 
(or more specifically, possibly admits a Killing tensor), 
and the absence of a constructive method to construct 
invariants. 

The two manifold approach to the problem has the 
great benefit that one can visually characterize the or- 
bits and identify the possibility of integrable behavior. 
It further allows one to illustrate the geometric meaning 
of a Killing tensor. One problem faced during calcula- 
tions is that the conformal factor J is very complicated, 
making the complete characterization of spacetimes for 
a given metric a formidable task and the characteriza- 
tion of all SAV spacetimes nearly impossible, using this 
approach. 

A large class of SAV spacetimes have orbits that ap- 
pear numerically to admit a fourth order invariant. This 
fact and the possibility of direct observational applica- 
tion if it does (Paper I of this series ) , has motivated 
a more in depth study (Paper IV Q). It turns out that in 
the context of the SAV spacetimes it may be possible to 
formulate an algebraic check that will determine whether 
a particular spacetime admits a higher order Killing ten- 
sor and thus quantify the relationship between the non- 
local metric distortion on a SAV spacetime described by 
the Weyl tensor and the dynamical behavior of particle 
motion within the spacetime. This formulation however 
requires the full power of the tetrad formalism and the 
solutiongeneration techniques, which I will review in Pa- 
per IV [a]. 
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APPENDIX A: SURFACES OF EXTREME 
ORBITAL CURVATURE 



2.2 



1.B - 



1.B - 



1.4 - 





2.4 2.5 2.6 27 2.8 2.9 



FIG. 7: Points of extreme orbital curvature in Schwarzschild 
E — 0.95, L = 3, i-i — 1. Black dots correspond to n = ±1 
surfaces. Blue have n = 0, 2 in equation (|A3|I . 

Points where an extremum of orbital curvature is 
reached can be computed by setting k(t) = 0. After 
some algebra one obtains 

k = [k, H] = ^y^j3/2 (°^(cos^ - sin^ 0) + 2b cos 9 sin d) , 

(Al) 

where a and b are functions of [^"^ , E,L, p, z) defined as 
follows: 

3 

a = -dpJdzJ - JdpzJ, 

b = \j{dpp d.,)J - \{{dpjf ~ [d.jf). (A2) 

As a result, curves in configuration space on which 
an extremum in curvature is reached [k — 0) can be 
parameterized by the phase angle 9c at the point where 

6lc = iarctan(^-^) ?i = ±0,l,2. (A3) 

The function 9c{p, E, L, jj?) has four possible solution 
surfaces with 9c G (— tt, tt] (depicted in figure [5]) This 
function can be thought of as a phase angle surface on 
which all points of extremal curvature must fall. 

For each surface a branch cut occurs if 6 = 0. Curves 
of extreme orbital curvature provide an accurate way of 
quantifying where, given a particular J, a particular orbit 
will be confined in configuration space. In effect what is 
constructed is a coordinate system ideally suited to the 
orbits. In FigurelHlthe extreme orbital curvature lines for 
several orbits for the J given in Figure [2] are computed 
and the four points of contact with the J = contour 
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Different branches of "S. Surfaces 




FIG. 8: The four branches on which solution points with ex- 
tremal curvature can lie 6c phase E = 0.95, L = 3, = —1. 



APPENDIX B: KILLING TENSORS IN 2 D 

The correspondence between the components of a 
phase space expansion Qn for calculating the invariant 
Q (equation I22p and polynomial in momenta and com- 
ponents of a Killing tensor on the two manifold is given. 

Consider the two metric gj^. = JSij on a two manifold 
admitting a totally symmetric Killing tensor - 
of order N. 

The Killing equations 

imply that 

Qt-T("--""W---P„„ (B2) 

remains constant along the geodesic. Making use of the 
definition of the momentum phase angle, equation (|16p . 
the invariant Qt can be rewritten as 




N 



FIG. 9: Lines of extreme orbital curvature for several orbits 
of a given J. Schwarzschild E = 0.95, Lz = 3, fj, = 1 



shown. Orbits on the left, that are not bound by a four 
pointed box, plunge through the horizon. 

Many quantities associated with the curvature and 
with extreme curvature surfaces can be most compactly 
expressed in complex notation. Introduce the complex 
variable = l/2(p + iz) and the complex function 
c = b + ia. Using equation (jA2p c can be expressed in 
terms of the the potential J as follows 



1=1 



(1 



(B3) 



where ai = 1 indicates the index p and aj = 2 the index 



z. To put this in the form of equation 



let P, 



{n,N) 



denote the j-th entry of the ?-th permutation of a list of 
a total of TV elements containing N — n entries equal to 1 
and n entries equal to -1. Let p{n) denote the number of 
permutations. As an example of the notation, let N = A. 
Then P(o,4)(l, j) = (1, 1, 1, 1) has one permutation, P(i,4) 
has 4 permutations, and P(2,4) has 6; these are explicitely 
listed below 

^(l,4)(l,j) = (-1,1,1,1) P{IA)%]) = (1,-1,1,1) 
^'(1,4)(3, j) = (1, 1, -1, 1) P(l,4)(4, j) = (1, 1, 1, -1) 
P{2A){l,j) = (-1,-1,1,1) -P(2,4)(2,j) = (-1,1,-1,1) 
^'(2,4)(3,j) = (-1,1,1,-1) P(2,4)(4,j) = (1,-1,-1,1) 
^(2.4)(5,j) = (1,-1,1,-1) P(2,4)(6,j) = (1,1,-1,-1) 

(B4) 

For clarity assume N is even. (The result also holds 
for odd N but more care has to be taken with the index 
ranges). The product of cosines can then be expressed 
as: 

N 



c = \jdaJ - I idcJf = -J"'da [-^] , (A4) 



and the extreme curvature conditions is 



1 

2^ 



N/2 P(^) (n 

E E e--expK:^/3,P(, 

l=-N/2 k=l \J = 1 



,N) 



(B5) 



(A5) 



Thus the correspondence between the even terms in 
the series for even N and the Killing tensor components 
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{21 



(B6) 



In the case where N = 21, the analytic function, men- 
tioned in equation ()25|) . is 



N 



= 2^T("i-"")i(^f=i(^~"^'). (B7) 



■N 2^ 



The exphcit expressions for the analytic function and 
Qq term of the lowest order case N = 2 are 



Qo = J{TPP + T''). 



(B8) 



For the = 4 case the expansion terms Qn expressed as 
a sum of the fourth order Killing tensor components are 



Manko-Novikov Spacetime 



R 



Vacuum 







Event 
Horizon 




Any set of mass multipole 
moments and spin can be 
specified. 

Horizon Area is 
explicitly determined 
by multipole moments 



Curvature Singularity on 

Horizon Equator 

(No singularities outside) 



FIG. 10: Properties of the Manko-Novikov Spacetime. 



— \{^PPPP ^ rpZZZZ _ QrjiZZpp _|_ ^^^rppppZ _ rppZZZy^ 
8 

Q2 = Y (TPPPP - T^^^^ - 2i{TPPP'' + TP"-'-)) , 



— (^J^PPPP _|_ 2T''-PP + J'zzzz-^ 



(B9) 



APPENDIX C: MANKO-NOVIKOV METRIC 

The metric of the Manko-Novikov spacetime V\ used to 
generate the plots in Figures [5l [6] and whose properties 
are sketched in Figure [TOl can be generated from the 
Ernst potential of the form 8 = e^"^ / where 



a;(l + ab) + iy{b — a) =F (1 — ia){l — ib), 



A^P = 0, 



(CI) 



and the coordinates x — coshp and y — cosz are called 
the Weyl coordinates. The metric functions that enter 
equation JT]) are 



The functions a and b obey a set of differential equations 
stated in [7|, and one example of a solution is given below. 

The solution for which the plots are made is a 
quadrupolar spacetime. Let r = {x^ + 2/^ — 1)^^^ and 
u — xy /r, and define the required Legendre polynomials 
to be Po{u) = 1, Pi{u) = u, P2{u) = -l/2 + (3/2)u2 and 
P^{u) = u (5u2 - 3) 12. Then 



In 



a 

—a 
b_ 
a 



In ( I ~2a2 

2a2 
P2 



i.x-y)[ — 
r 



Pj. 

.3 



{x + y) 



Po 



tp — a2 

I- 

1 , 



- y 
' p.? 



p2 
^2 



— +ln - 



(C4) 



k'ix'-m-y'), e 



where 



27 , Mx'' - y^) 

(a;2- l)(l-a2)2^ 



A 1 - a2 ' 
(C2) 



A={x^~ 1)(1 -I- a6)2 - (1 - y^){b - a)2, 
B=[x + l + {x- l)a6]2 + [(1 + y)a + (1 - y)bf , 
C^{x^- 1)(1 + ab)[b -a-y{a + b)] 

+ {l-y^){b-a)[l + ab + xil-ab)]. (C3) 



Note that ip and 7 are members of the Weyl class of static 
metrics, a is a parameter that scales the spin and 012 is 
the quadrupole moment. The Geroch-Hanson multipole 
moments for this metric can be found in Q along with 
a more general solution parameterized by arbitrary mass 
multipole moments. In the event that a = and a2 = 
the metric reduces to the Schwarzschild metric with met- 
ric functions 



i?2 



1 



X+1^ 

(x2-l)(l-y2), 



.27 



(jj 



X 



X'^ - 

0, k 



(C5) 
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APPENDIX D: ZIPOY-VOORHEES METRIC 

The Zipoy-Voorhees metric [13] is a static space- 
time with metric functions 

i?2 (a;2 - 1)(1 -2/2)^ ^ = 0, k=l. (Dl) 



All numerical experiments performed in this metric thus 
far appear to admit integrable orbits similar to that por- 
trayed in Figure m 
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